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With a sample of 232 xlO 6 T(AS) -» BB events collected with the BABAR detector, we study 
the decay B + — > ppK + excluding charmonium decays to pp. We measure a branching fraction 
B(B + — > ppA" + ) = (6.7±0.5±0.4) x 10 -6 . An enhancement at low pp mass is observed and the Dalitz 
plot asymmetry suggests dominance of the penguin amplitude in this B decay. We search for a 
pentaquark candidate decaying into pK + in the mass range 1.43 to 2.00GeV/c 2 and set 

limits on B{B+ -> 6*++p) x B{B*++ -> pK+) at the 10" 7 level. 

PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 



This paper describes a measurement of the branch- 
ing fraction of the baryonic three-body decay B + — > 
ppK + (excluding charmonium decays to pp) and a 
study of its resonant substructure. An earlier measure- 
ment of the branching fraction for this channel gave 
(5.7tg;g± 0.7) x 1CT 6 . This channel is interesting for the 
dynamical information in the distribution of the three 
final-state particles and for the possible presence of ex- 
otic intermediate states. We also isolate decays 
B + — > X C c.K + , where X cS = r\ c and J/ip decaying to pp, 
and measure the width of the rj c . 

An important feature of this decay is an enhancement 
at low pp masses reported in Ref. , similar to those that 
have been observed in several other baryonic decays of 
B and J/ip @. This could be a feature of a quasi- two- 
body decay in which the pp system is produced through 
an intermediate gluonic resonance Q (Fig-QJc)). It could 
also be that the decay is a pure three-body process and 
that the enhancement results from the short-range corre- 
lations between p and p in the fragmentation chain 0] > @ • 
Rosner suggested that it is possible to distinguish 
the fragmentation or gluonic resonance mechanisms by 
studying the distribution of events in the Dalitz plot. 

The main Feynman diagrams for this decay are pre- 
sented in Fig. n] The leading diagrams are a penguin 
diagram and a doubly Cabibbo-Kobayashi-Maskawa- 
suppressed tree diagram shown in Fig. ^a,b). There 
is also an Okubo-Zweig-Iizuka-suppressed penguin dia- 
gram shown in Fig. ^c), where the pp pair is created 
through a pair of gluons (or a gluonic resonance). There 
are four additional color-suppressed diagrams ; 8]: two 
tree diagrams with an internal PF + -emission and a W + - 
annihilation and two penguin diagrams with an internal 
gluon-emission that are expected to be small. If the pp 
system is produced independently of the K + through a 
tree diagram with an external IF + -emission (Fig. E^b)) 
or a penguin with an external gluon-emission (Fig. 
i.e. the pp quark lines are not associated with the s or 
u quarks in the K + , then the distributions m pK + and 
rripK+ should be identical. If the u quark in the K + is 
associated with a u quark in a p (Fig. Ufa)), larger val- 
ues of m pK + are favored over those of rrip K + @ . Thus a 
study of the Dalitz plot provides insight not only into the 




FIG. 1: The main Feynman diagrams for the non-resonant 
B + — » ppK + decay: (a) leading penguin diagram, (b) lead- 
ing tree diagram (external W^-emission), (c) Okubo-Zweig- 
Iizuka-suppressed penguin diagram. 

dominant mechanism of this decay but also into whether 
the penguin or the tree amplitude is dominant. 

This paper is organized as follows: first we describe 
the event selection and the branching-fraction measure- 
ment. Then we describe the pp mass spectrum and the 
measurement of the r) c width. We examine the Dalitz 
plot for an asymmetry between the distributions in m p K+ 
and rripK+ . In the final section we describe searches for 
B + — > pA(1520) — > p(pK + ) decay and for the hypothe- 
sized I = 1, J3 = 1 pentaquark state 0* ++ (a member of 
the baryon 27-plet with quark content uuuds) in the de- 
cay B + — > Q* ++ p —> (pK + )p. The 9* ++ mass has been 
predicted 10] to lie in the region 1.43 — 1.70GeV/c . 

The analysis uses 232 xlO 6 Y{AS) BB decays col- 
lected with the BABAR detector [ill at tnc PEP-II e+e" 
storage ring. Charged tracks are measured by a five-layer 
silicon vertex tracker (SVT) and a 40-layer drift-chamber 
(DCH) in a 1.5-T solenoidal magnetic field. A Cherenkov 
radiation detector (DIRC) is used for charged-particle 
identification. The CsI(Tl) electromagnetic calorimeter 
detects photon and electron showers. To identify kaons 
and protons we use dE/dx measurements in the SVT 
and DCH, and the pattern of Cherenkov photons in 
the DIRC. The proton efficiency is 93% with 9% kaon 
misidentification probability. The kaon efficiency is 87% 
with 2% pion misidentification probability. 
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We use the kinematic constraints of £?-meson pair- 
production at the T(4S) to identify the B + — > ppK + 
signal. Two independent variables are calculated for each 
ppK+ candidate: m ES = [(£ 2 m /2+p -p B ) 2 / Efi-pl] 1 / 2 
and AE = Eg — E cm /2, where E cm is the total center- 
of-mass energy, the subscripts and B refer to the initial 
T(4S) and to the B candidate, respectively, and the as- 
terisk denotes the T(45) frame. The resolutions on AE 
and mES are about 17MeV and 2.6MeV/c 2 , respectively. 

Backgrounds arise primarily from random combina- 
tions in continuum events (e + e~ — * qq 1 where q = 
u,d,s,c). These events are collimated along the origi- 
nal quark directions and can be distinguished from more 
spherical BB events with a Fisher discriminant (J 7 ) ^3] j 
a linear combination of four event-shape variables. The 
four variables are cos9^ hr , the angle between the thrust 
axis of the reconstructed B and the beam axis; cos6^ om , 
the angle between the momentum of the reconstructed 
B and the beam axis; and the zeroth and second Legen- 
dre polynomial momentum moments, Lq = ^\ |p*| and 
L 2 = E 4 |P 2 *l[(3cos 2 6* thrB i - l)/2], where p* are the mo- 
menta of the tracks and neutral clusters not associated 
with the B candidate and Q* thTB . is the angle between p* 
and the thrust axis of the B candidate. The event selec- 
tion is optimized assuming the previously measured value 
of the branching fraction |2j to maximize s/ys + b, where 
s and b are the expected number of signal and back- 
ground events, respectively. The event selection retains 
66% of signal events while removing 93% of continuum 
background. 

The resulting distribution of events in the mES~AE 
plane is shown in Fig. [21 A clear signal is ob- 
served at the B mass and AE=0. Potential back- 
grounds are studied with Monte Carlo (MC) simula- 
tion [l4|. The combinatorial background is expected 
to come predominantly (89%) from continuum events. 
Background events in the signal region arise mostly 
from B + — > X c5 (pp)K + , where X c£ = ri c , J/ip, tp', Xc0,l,2 
(the charmonium background), while non-charmonium 
B backgrounds are expected to be negligible. The 
signal and sideband regions are defined to be "wide" 
(5.27 < m ES <5.29 GeV/c 2 and 5.20<m B5 <5.26 GeV/c 2 , 
|A_E|<5GMeV) for the charmonium background stud- 
ies and "narrow" (5.276<m£s<5.286 GeV/c 2 and 
5.20<m£;5<5.26GeV/c 2 , | A£|<29 MeV) for the Dalitz 
plot study. 

To extract the ppK + signal yield, we fit the 
AE distributions for candidates that lie in the 
5.27<m£;s<5.29 GeV/c 2 region separately in nine bins 
of m p p. The size of the bins is shown in Fig. [3J We use 
a linear function for the background and a double Gaus- 
sian distribution for the signal. The widths and means 
of the Gaussian distributions and their relative areas are 
fixed to values obtained from MC simulation, which is 
also used to calculate the detection efficiency {s mpP ) in 
each m„p bin. Across the allowed kinematic region, e m 5 
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FIG. 2: Distribution of AE versus tues for the ppK + candi- 
dates in data. The solid (dashed) lines define the wide (nar- 
row) signal and sideband regions. 
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FIG. 3: The m p p distribution for data in the wide signal 
(points) and sideband (shade) regions. The sideband his- 
togram is scaled to the expected number of the combinatorial 
background events in the signal region. 

declines smoothly from 30% at threshold to 24% at the 
highest mass. The AE fits for m p p below 2.85GeV/c 2 
yield 343l 2 6 signal events. From the known number of 
charged B mesons in the sample, the branching fraction 
for m p p below the rj c mass is measured to be B{B + — > 
ppK+- m P p<2.85 GeV/c 2 ) = (5.3 ± 0.4 ± 0.3) x I0~ 6 . 

An estimate of the number of charmonium events 
in the m P p>2.85 GeV/c 2 region is required to deter- 
mine the total non-charmonium branching fraction. To 
minimize the systematic error on that quantity, we 
fit the m p p spectrum for the number of the non- 
charmonium events in the primary "charmonium" region 
(2.85<m P p<3.15 GeV/c 2 ). To improve the pp mass res- 
olution in the m pp fit, we perform a kinematic fit fixing 
the mass and energy of each B candidate in the wide 
signal and sideband regions to their known values. The 
m pp distribution is shown in Fig. where prominent 
signals for the r\ c and J/if) decaying into pp are visible. 
The region used in the m pp fit, 2.4<m P p<3.4 GeV/c 2 , 
is chosen wider than the "charmonium" region of inter- 
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est (2.85<m P p<3.15 GeV/c 2 ), shown in Fig. inset) , to 
improve the statistical uncertainties on the ppK + signal 
and combinatorial background yield. The r\ c peak is de- 
scribed by a convolution of a Breit-Wigner distribution 
and a Gaussian distribution, and the J/ip peak by a sum 
of two Gaussian distributions with a common mean. The 
shapes are obtained from MC simulation. The width of 
the broader J/ip Gaussian distribution and ratio of areas 
of the two J I ip Gaussian distributions are constrained in 
the fit to their MC values. A common width is used for 
the narrow Gaussian distributions for J/ip and r\ c and is a 
free parameter in the fit. The ppK + signal and combina- 
torial background distributions are modeled by a linear 
function of m p p. The inset of Fig.|2lshows this fit, which 
results in 1147^ rj c events and 1371^ ■VV' events. Cor- 
recting for the detection efficiency of (26.9 ± 0.2)%, we 
find 23(5+ -> rj c K+)xB(r] c ^ pp) = (1.8±^±0.2)x 1CT 6 
&ndB(B+ -> J/ipK+)xB(J/ip ->■ pp) = (2.2±0.2±0.1)x 
10~ 6 in agreement with the accepted values The fit 
yields a total rj c width of r(?7n)=25±5±3 McV/c 2 consis- 



tent with the current values 13] and a mass resolution of 
5.7 ± 0.4MeV/c 2 in agreement with MC expectations. 

The m p p fit yields 88±6 ppK + signal and combina- 
torial background events in the charmonium region (see 
Fig. In this region, the latter contribution is esti- 
mated from the AE fit to be 53±5 events, resulting in 
a non-charmonium ppK + signal of 35±8 events. The 
contribution from higher-mass charmonium modes is es- 
timated to be 24±5 events from the accepted 0] val- 
ues for their branching fractions. Adding the ppK + sig- 
nal yield obtained from the AE fits outside the "char- 
monium" region with non-charmonium ppK + signal in 
the "charmonium" region, and subtracting the contribu- 
tion of the higher mass charmonium modes results in 
a total non-charmonium signal yield of 433±33 events. 
Correcting the signal yield for efficiency in each of the 
m p p bins and normalizing to the number of B + mesons 
in the data sample results in a total branching fraction 
of B(B+ -> ppK+) = (6.7 ± 0.5 ± 0.4) x 10~ 6 with 
charmonium decays to pp excluded. Figure 0] shows the 
background-subtracted and efficiency-corrected pp mass 
spectrum and the expectation for a three-body phase- 
space decay. The existence of a low-mass enhancement 
in the pp mass as previously observed by Belle 2j is con- 
firmed. 

The charge asymmetry is defined as 
A ch =(N B —N B +)/(N B -+N B +), where N B ± is the 
number of B^—^ppK^ events. We use the same fitting 
procedure as for the branching fraction measurement, 
and find A cfc =-0.16±g;^ ± 0.04 for m pp <2.85 GeV/c 2 . 

For the remainder of this paper to increase the sig- 
nal purity, only events in the narrow signal and ities- 
sideband regions are considered. After selecting the B 
candidates, we perform a kinematic fit for each B candi- 
date, fixing its mass and energy to their known values. 

We study the dynamics of the three-body decay by 
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FIG. 4: Efficiency-corrected yield of B + — » ppK + events as 
a function of m p p in data (points) and in three-body phase- 
space signal MC (histogram). Errors shown are statistical 
only. 
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FIG. 5: Dalitz plot of data in the narrow signal region, and 
sideband region (inset). These distributions are not efficiency- 
corrected. The lines show the positions of the prominent char- 
monium backgrounds, from left to right ip', Xc2,i,o, J/ip, rjc- 
The sideband contains about eight times more combinatorial 
events than are expected in the signal region. 

constructing signal and sideband Dalitz plots (Fig. [SJ. 
There are 780 (1661) events in the signal (sideband) re- 
gion. The sideband contains about eight times more com- 
binatorial events than the signal region. The Dalitz plot 
for the signal region shows the threshold enhancement in 
the pp mass spectrum, as well as clear diagonal bands 
corresponding to r/ c , J/ip and ip 1 decays. 

To study the m p x+ and rn p K+ asymmetry, we divide 
the Dalitz plot along the m p K+ = m pK + line (dashed 
line in Fig. [SJ and each of the two halves is projected 
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FIG. 6: Efficiency-corrected distributions in the narrow sig- 
nal (points) and rescaled sideband (shade) regions: (a) m pK + 
(for m pK + > m pK + ) , (b) m pK + (for m pK + < m pK + ) , and (c) 
difference between (a) and (b). Errors shown are statistical 
only. 



onto the nearer axis. The corresponding distributions 
for the events in signal and rescaled sideband regions are 
shown in Fig. Ef a 5 b). No asymmetry is expected from 
variations in s mpP which is charge-symmetric and slowly 
varying with m p p, nor from the small combinatorial back- 
ground shown in Fig.|^a,b). The asymmetry appears as 
a broad enhancement peaking at about 4 GeV in the pK + 
combinations (Fig. Etc))- This could be an indication of 
a correlation between quarks in p and K + if the B de- 
cay proceeds through a penguin diagram (Fig.^a)). No 
quantitative theoretical description of this correlation is 
available at the moment. 

The two-body decay B + — > p4(1520) could also be 
present in the ppK + signal sample. The efficiency for 
detection of this channel is determined in dedicated MC 
simulation to be (4.7 ± 0.1)%, including ,8(4(1520) -> 
pK~) 13]. The mpK+ spectrum, shown in Fig. Ufa), is 
fit with an ARGUS function [l5j for the background and 
a Breit-Wigner convolved with a double Gaussian (with a 
common mean) for the 4(1520) signal shape. The mass 
resolutions and the ratio of areas of the Gaussians are 
fixed to the values obtained from MC simulation, while 
we fix the mean and the natural width to established 
values 0] ; the endpoint of the ARGUS function is fixed 
to the sum of the proton and kaon masses. An unbinned 
maximum likelihood fit (Fig. da)) results in an upper 
limit (U.L.) on B{B+ -> p4(1520)) of 1.5 x 10~ 6 at 90% 
C.L. (including a systematic error of 16%). 

The search for a ligh t pentaquark candidate 

(m e *++<2 GeV/c 2 ) [lfj proceeds as follows. From 
B + — > ppK + three-body phase-space MC as well 



TABLE I: Systematic uncertainties in percent on the branch- 
ing fraction measurements and in the values of uncertain- 
ties for the symmetry measurements. Values for m p p below 
2.85GeV/c 2 are given in brackets. 



Type 


ppK + 




pA(1520) 


p6* ++ 


A ch 


B-counting 


1.1(1.1) 


1.1 


1.1 


1.1 




Tracking/PID 


3.8(3.8) 


3.4 


4.2 


4.2 


0.02 


MC Statistics 


2.1(2.4) 


0.7 


1.0 


0.5 


0.03 


B.F. Errors 


0.9(-) 




2.2 






Selection 


0.2(-) 


0.4 


3.9 


3.9 




AS/Mass Fits 


3.6(2.4) 


8.9 


14.3 




0.01 


Total 


5.8(5.2) 


13.5 


15.6 


6.1 


0.03 




, & 4.5 



H 3.5 
C 



> 



(b) 



1.4 1.45 1.5 1.55 1.6 1.65 1.7 1.75 



m_ ,(GeV/c ) 

pK 



m pK .(GeV/c) 



FIG. 7: (a) The m pK + distribution for data events in /1(1520) 
mass region; (b) The m pK + distributions for data events 
in the signal region outside (dashed) or inside (solid) the 
2.85<m P p<3.15 GeV/c 2 region. Distributions are not effi- 
ciency corrected. 

as five dedicated signal MC samples with 77i e *++ = 
1.5, 1.6, 1.7, 1.8, 1.9GeV/c 2 , we find the mass reso- 
lution {(J p K+) t° vary from 1.0 to 3.0MeV/c 2 for 
1.43<m p / f +<2.00 GeV/c 2 , and the average reconstruc- 
tion efficiency to be (20.5 ± 0.1)% in this mass region. 
The events with m p p in the charmonium region are ve- 
toed. The pK + mass distribution of the remaining 
events is shown in Fig. EJb)- A Bayesian approach is 
used to calculate the U.L. at 90% C.L. as a function of 
m pK+i assuming Poisson-distributed events in the ab- 
sence of background. Each limit is increased by the 
total systematic error of 6%. The U.L. for B{B + — > 
e*++p) x B{Q*++ -> pK+) is measured to be 0.5 x 10~ 7 
for 1.43 < m(e*++) < 1.50GeV/c 2 , < 0.9 x 10~ 7 for 
1.50 < 77i(e*++) < 1.72GeV/c 2 , and < 1.2 x 10~ 7 for 
1.72 < m(9*++) < 2.00GeV/c 2 . 

The systematic uncertainties for each analysis are sum- 
marized in Table ||| The T(4S) is assumed to decay 
equally to B°B° and B + B~ mesons. Incomplete knowl- 
edge of the luminosity and cross-section leads to a 1.1% 
uncertainty. Charged-tracking and particle-identification 
(PID) studies in the data lead to small corrections ap- 
plied to each track in these simulations. Limitation of 
statistics and purity in these data-MC comparisons lead 
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to residual tracking/PID uncertainties. A large control 
sample of B + —> J/ip(e + e~)K + is separately studied in 
data and MC simulations to understand the residual er- 
rors from the event-shape, AE, and titles cuts. Lim- 
itation of MC statistics employed in each analysis con- 
tributes to a small uncertainty. Branching fraction uncer- 
tainties (B.F. Errors) on B{B+ -> XK+) x B{X -> 
pp), where X = Xc[o,i,2] > ^ an( i 8(^1(1520) — > pK~) af- 
fect the total ppK + and the pA branching fraction mea- 
surements, respectively. Where the MC values are used 
to fix signal shape parameters in a fit, the parameters are 
varied within their uncertainties and the data are refit to 
propagate this uncertainty. In a similar fashion, different 
ranges and background functions are employed to estab- 
lish the uncertainty on the mass spectra fits (resulting, 
for example, in the T(?7 C ) uncertainty of 3MeV). 

In summary, with 210 fb -1 of data, we isolate the 
B + — > ppK + final state, and measure its non- 
charmonium branching fraction to be (5.3 ± 0.4 ± 0.3) x 
10~ 6 for m pP below 2.85GeV/c 2 and (6.7 ± 0.5 ± 
0.4) x 10~ 6 for the whole uipp range. We measure 
A c fc=-0.16lg;^ ± 0.04 for m pP below 2.85GeV/c 2 . The 
existence of a low-mass enhancement of the pp pair is 
confirmed. The asymmetry between pK + and pK + fi- 
nal states in the Dalitz plot is demonstrated, providing 
evidence supporting the dominance of the penguin am- 
plitude in this B decay. We measure the total width of T] c 
to be 25^5±3 MeV/c 2 . An upper limit of the decay rate 
to pA(1520) is set at 1.5 x 10~ 6 . No evidence is found 
for the pentaquark candidate 9* ++ in the mass range 
1.43 to 2.0GeV/c 2 , decaying into pK + , and branching 
fraction limits are established at the 10 -7 level. 
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